arXiv:1506.01754vl [physics.optics] 5 Jun 2015 


Optical forces in nanorod metamaterial 


Andrey A. Bogdanov 1,2 Q Alexander S. Shalin 1 , and Pavel Ginzburg 3 
1 ITMO University, 197101 St. Petersburg, Russian Federation 
2 Ioffe Institute, 194 021 St. Petersburg, Russian Federation 
3 School of Electrical Engineering, Tel Aviv University, 69978 Tel Aviv, Israel 

(Dated: June 8, 2015) 

Optomechanical manipulation of micro and nano-scale objects with laser beams finds use in a 
large span of multidisciplinary applications. Auxiliary nanostructuring could substantially improve 
performances of classical optical tweezers by means of spatial localization of objects and inten¬ 
sity required for trapping. Here we investigate a three-dimensional nanorod metamaterial platform, 
serving as an auxiliary tool for the optical manipulation, able to support and control near-field inter¬ 
actions and generate both steep and flat optical potential profiles. It was shown that the ’topological 
transition’ from the elliptic to hyperbolic dispersion regime of the metamaterial, usually having a 
significant impact on various light-matter interaction processes, does not strongly affect the distri¬ 
bution of optical forces in the metamaterial. This effect is explained by the predominant near-fields 
contributions of the nanostructure to optomechanical interactions. Semi-analytical model, approxi¬ 
mating the finite size nanoparticle by a point dipole and neglecting the mutual re-scattering between 
the particle and nanorod array, was found to be in a good agreement with full-wave numerical simu¬ 
lation. In-plane (perpendicular to the rods) trapping regime, saddle equilibrium points and optical 
puling forces (directed to the excitation light source along the rods), acting on a particle situated 
inside or at the nearby the metamaterial, were found. The auxiliary metamaterials, employed for 
optical manipulation, provide additional degrees of freedom in flexible nano-mechanical control and 
could be employed in various cross-disciplinary applications. 


I. INTRODUCTION 

The ability to control mechanical motion of micro- and 
nano-scale particles with focused laser beams is an essen¬ 
tial tool, being a paramount for a wide range of applica¬ 
tions, related to bio-physics, micro-fluidics, optomechani¬ 
cal devices and more! 1 ^ Being first proposed and demon¬ 
strated by A. AshkinP the classical optical tweezers are 
nowadays a rapidly developing area of fundamental and 
applied research. 

One of the promising and already conceptually proven 
approaches for improving performances of the optical ma¬ 
nipulation schemes is to employ various auxiliary nanos¬ 
tructures, especially plasmonic ones? 1 The key idea of 
the plasmonic tweezers is to utilize strong light-matter 
interaction between nanostructured metals and focused 
laser beams. Noble metals, having a negative permittiv¬ 
ity in the optical and infrared spectral ranges, support 
localized plasmon resonances enabling enhancement and 
control of near-fields at their vicinity. 6 7 In particular, 
the creation of strong intensity gradients is beneficial for 
obtaining substantial optical forces, which is important, 
for example, for achieving molecular manipulation.® 

Plasmonic nanostructures with subwavelength light 
concentration could be employed for obtaining new op¬ 
tomechanical effects, i.e. accelerating nanoparticles in an 
arbitrary direction (in relation to the light propagation 
direction), 9 or for creating nano-modulators of plasmonic 
signals 

Array s of antennas and their integrations in photonic 
circuitryp^J employed as auxiliary tools for optical trap¬ 
ping, were shown to outperform classical schemes (fo¬ 
cused lasers in homogeneous media, e.g. liquid solu¬ 


tions) both in terms of spatial localization and optical 
power required per trapped particle.Antenna arrays were 
further extended for multifunctional platforms, enabling 
trapping, stacking, and sorting.^ However, isolated plas¬ 
monic structures create limited number of hot spots (lo¬ 
cal enhancement of intensity) and are usually restricted 
to two-dimensional geometries. These constrains set sig¬ 
nificant limitations on the flexibility of optical manipula¬ 
tion by reducing potential degrees of freedom, available 
for optomechanical control. On the other hand, three- 
dimensional artificially created nanostructures or meta¬ 
materials (see, e.g., Ref. fl3l) could provide additional ben¬ 
efits and flexibility by configuring near-field interaction 
in large volumes. 

Hyperbolic metamaterial^^ are one class of artificially 
created electromagnetic structures, capable to enhance 
efficiencies of various light-matter interaction processes 
owning to the unusual dispersion of eigenmodes, sup¬ 
ported by the structure - namely its hyperbolic disper¬ 
sion. Among various possible designs of this type of 
metamaterials it is worth m ention ing composites made 
of vertically aligned nanorods JI5UGJ periodic metal dielec¬ 
tric layersp^and semiconductor quantum structures! 18 * 19 ! 
While the far-field interactions of waves with hyperbolic 
composites were proven to be well characterized in terms 
of the effective medium approximation, this description 
could be questioned if near-field mediated processes are 
involved [see, e.g., Ref. EDI- 

The general criterion for estimating impact of near¬ 
field contributions to an interaction is based on compari¬ 
son of k-vector spectra with reciprocal vector of the meta¬ 
material lattice. For example, scattering from objects 
within hyperbolic metamaterials involves consideration 
of the near-field effects [see, e.g., Ref. EQ. 
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FIG. 1. (color online) Schematics of the nanorod metama¬ 
terial with a spherical nanoparticle inside. Radius of the 
nanoparticle is R = 5 nm. Height and radius of the nanorods 
are H — 350 nm and r — 15 nm, respectively, the period is 
h — 60 nm. 
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FIG. 2. (color online) Frequency dependence of the effec¬ 
tive tensor components (real parts) £± and £\\ of the nanorod 
metamaterial. Dashed line shows the transition between el¬ 
liptic and hyperbolic dispersion regimes of the metamaterial 
- epsilon near zero (ENZ) point. Insets show characteristic 
shapes of iso-frequency surfaces, corresponding to the disper¬ 
sion regimes. 


Analysis of near- and far-field contributions to opti¬ 
cal force, acting on objects embedded in the nanorod 
metamaterial is the central topic of the manuscript. In 
particular, optical forces, acting on nano-sized spheri¬ 
cal particle embedded inside the metamaterial assembly, 
are investigated both numerically and by using a semi- 
analytical approach, considering the finite size nanoparti¬ 
cle as a point dipole and neglecting re-scattering between 
the particle and nanorod array. The impact of the finite 
structure of the metamaterial unit cell and the relative 
arrangement of the particle in respect to it was analyzed 
as a function of the system’s geometry and frequency of 
incident illumination. A semi-analytical approach based 
on dipole near-field interaction is developed and shown 
to be in a good agreement with results of the full-wave 
numerical analysis. The interplay between near- and far- 
field effects in the context of effective medium approxi¬ 
mation is discussed. 

The manuscript is organized as follows. In Sec. [II] the 
nanorod metamaterial platform is introduced and its op¬ 
tical properties are analyzed. In Sec. m the distribu¬ 
tion of optical forces is investigated numerically. A semi- 
analytical model is introduced and compared with the 
full-wave numerical simulation. Finally, Sec. |TV| summa¬ 
rizes the major results. 


II. NANOROD METAMATERIAL: FAR-FIELD 
CHARACTERISTICS 

The geometry under investigation is schematically rep¬ 
resented in Fig. [l]- it shows an array of vertically aligned 


gold nanorods and a gold nanoparticle placed inside it. 
Material parameters of the constitutive elements were 
taken from widely used sources.^ The parameters of the 
structure are indicated in the figure caption. While the 
nanorods in this model are situated in vacuum, sub¬ 
strate effects and host material filling the space be¬ 
tween the rods could be taken into account straight¬ 
forwardly. Similar structures have already found use 
in various multidisciplinary applications, among them 
bio-sensingP^ enhancement of nonlinearitiesacous¬ 
tic waves detection,^ thin optical elemental and oth¬ 
ers. The key properties of this auxiliary nanostructure 
leading to enhanced performance are large surface area 
and unusual collective optical response of the system, 
enabling control over both far- and near-field interac¬ 
tions. Hence, investigation of optical forces, mediated 
by nanorod metamaterials, has a profound potential in¬ 
terest. 

Far-field interactions between electromagnetic waves 
and metamaterials, under certain circumstances, can be 
described within the effective medium approximation. 
The main idea of this homogenization procedure is to 
average the electromagnetic field over a unit cell of the 
structure. Therefore, the field inside the structure is as¬ 
sumed to be uniform. In the context of optical forces, 
as it will appear in the next section, the non-uniformities 
play a major role and, in fact, predefine the spatial struc¬ 
ture of optical potentials. Recently, a phenomenological 
approach taking into account the finite size of the meta¬ 
material unit cell was proposed.^ Inclusion of a depo¬ 
larization volume around optically manipulated particles 
enabled investigations of far-field contributions to optical 
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forces. However, near-field interactions, being strongly 
dependent on a specific metamaterial design, were not 
included explicitly. 

One of the key properties of hyperbolic metamaterials, 
making them attractive for electromagnetic applications, 
is their unique ability to support an unusual regime of 
dispersion caused by having permittivity tensor compo¬ 
nents of opposite signs (e±£\\ < 0). An immediate im¬ 
plication of this hyperbolic dispersion regime is the high 
density of ph otonic states, available for both emission 

and scattering 

The effective permittivity tensor of nanorod metama¬ 
terial is given by: 

fe± 0 0\ 

e= 0 cjl 0 . (1) 

\0 0 V 

Here and ey are effective permittivities perpendicular 
and along the wires, respectively. 

Dispersion of the tensor components for the structure 
under consideration was calculated with the approach de¬ 
veloped in Ref. [29] The transition between elliptic and 
hyperbolic dispersion regimes occurs at the wavelength 
around 523 nm (see Fig.j2|. The transition point is called 
epsilon-near-zero (ENZ) regime, as the real part of the 
permittivity along the rods is vanishing, if the spatial 
dispersion effects are ignored. The high density of pho¬ 
tonic states as well as the strong scattering emerges in the 
hyperbolic and ENZ regimes. The wavelength of the ex¬ 
ternal illumination, exploited for optical manipulation in 
the subsequent investigations, is chosen around this ENZ 
point in order to distinguish between various dispersion 
regimes and their impact on optical forces. 

III. OPTICAL FORCE DISTRIBUTION 
A. Numerical model 

The distribution of the optical forces, acting on the 
gold nanoparticle placed inside the wire medium is ana¬ 
lyzed hereafter. Both the hyperbolic and elliptic disper¬ 
sion regimes of the bulk metamaterial and their impacts 
on optical forces are compared and discussed. 

The incident wave is chosen to be linearly polarized 
along the y -axis and it propagates along the z- axis, as 
it is shown in Fig. [T| — standard normal incidence sce¬ 
nario is considered. Full 3D numerical analysis, based on 
finite elements method,^ is performed in order to calcu¬ 
late self-consistent electromagnetic fields in the system. 
Consequently, optical forces acting on the nanoparticle 
are calculated by integrating the Maxwell’s stress ten¬ 
sor components over an imaginary spherical surface sur¬ 
rounding the nanoparticle. 

The presence of a single nanoparticle breaks the inher¬ 
ent translation symmetry of the initial metamaterial ge¬ 
ometry. In order to overcome the computation complex¬ 


ity of large systems modeling, Floquet periodical bound¬ 
ary conditions were imposed on finite size geometries. 
This type of model corresponds to a periodic system with 
variable unit cell, which consists of a square array of 
nanorods and the nanoparticle. If the electromagnetic 
coupling between the particles in adjacent cells is minor, 
this type of analysis recovers the behavior of the infinite 
system with a single particle. 

The numerical procedure is as follows: the number 
of rods in the unit cell is increased gradually and the 
convergence of a certain quantity (optical forces in our 
case) is checked. Recently, a similar approach was ap¬ 
plied in studies of the Purcell effect in nanorocPH and 
wir (P metamaterials. Square unit cells containing 4, 9, 
and 16 nanorods were considered and the convergence of 
optical force values at different points of the metamate¬ 
rial volume was checked. A unit cell of 4 nanorods (the 
smallest one) was shown to predict the behavior of an 
infinite array within the accuracy of several percent. All 
the subsequent results were obtained for this size of the 
unit cell. The direct consequence of this calculation is 
that (i) only nearest neighbor rods define the value of 
optical force and (ii) nanoparticles in different unit cells 
almostly do not interact with each other. 

It should be noted, however, that the collective macro¬ 
scopic behavior of the array is taken into account by im¬ 
posing periodical Floquet boundary conditions. 


B. Lateral force component 

All the subsequent calculations were done for a par¬ 
ticle of 10 nm in diameter. The optical force F, in 
the most general case, has three non-zero components 
(F x , F y , F z ). The lateral force = (F x ,F y , 0) will be 
analyzed first. Values of optical forces are normalized to 
the intensity of the incident wave and volume of the par¬ 
ticle in order to perform direct comparisons with other 
optical manipulation schemes. 

The resulting normalized forces at the cut-plane 
z = 10 nm, calculated for wavelengths A = 450 and 
600 nm, corresponding to the elliptic and hyperbolic 
dispersion regimes respectively, are shown in Fig. |3ja). 
Both cases are shown on the same panel, but with dif¬ 
ferent scales of the color bar, as the impact of the exci¬ 
tation wavelength on the spatial distribution of the op¬ 
tical forces is minor. Different color bars indicate the 
differences in the absolute values. Forces maps at var¬ 
ious cut-planes (with different ^-coordinate) show qual¬ 
itatively similar behavior. It should be noted, that the 
values and directions of forces attributed to the geomet¬ 
rical center of the particle, hence certain regions (dark 
blue shells around nanowires with thickness equal to the 
radius of nanoparticle) on Fig. |3ja) and (c) are blank, as 
this center cannot approach the boundaries of the rods. 

It can be seen from Fig. |3|a) that the force distribu¬ 
tion has saddle points at the center of the unit cell and 
at its edges between the rods. These places correspond 
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FIG. 3. (color online) (a) Lateral optical force distribution at the cut-plane z — 10 nm (see Fig. [lj. Magnitude of the 
lateral force F | = (F x ,F y , 0) is shown with a color scale and the direction with arrows, (b) Distribution of the electric field 
magnitude in the cut-plane z = 10 nm. Magnitude of the electric field is shown with the color scale. The arrow lines are 
electric field lines. Panel (c) shows the distribution of the force calculated with a semi-analytical dipolar approach [see Eq. ©]. 
i.e. where the perturbation of the field by the particle is neglected. Panel (d) shows the numerical simulation of the electric 
field magnitude distribution without the particle. Panel (a) should be compared with (c), while (b) with (d). Dark blue shells 
around nanorods on panels (a), (c) have the width equal to particle’s radius. Optical forces are not calculated at those areas, 
as the nanoparticle’s center cannot approach the nanorods that close. Upper and lower scales of the color bars correspond to 
the hyperbolic (A = 600 nm) and elliptic (A = 450 nm) dispersion regimes of the metamaterial, respectively. Electric field 
amplitude of the incident wave is 1 V/m. 


to the the saddle points of electromagnetic field magni¬ 
tude distribution [Fig. J3|d)], and, consequently, to the 
unstable equilibrium positions of the nanoparticle. Some 
peculiarities in the optical force distribution appear on 
cut-planes near the edges of the nanorods (z = 350 nm 
and z = 0 nm), but they are attributed to the longitudi¬ 
nal (z-component) force component and will be discussed 
in Sec. HITO1 

The similarity of the spatial distribution of the forces 


at hyperbolic and elliptic dispersion regimes results from 
the dominating near-field coupling between the nanorods 
and the particle. Figure |3][b) shows the magnitude dis¬ 
tribution of total electric field |E| = (E% + Ey + E J) 1 / 2 , 
while the arrows show its direction. One can see that the 
field map is formed by electrical dipoles induced on the 
rods and the particle by the incident wave. Orientations 
of the dipoles nearly coincide with the polarization of the 
incident wave. Minor deviations from the above descrip- 
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FIG. 4. (a) Distribution of: (a) Total electric field magnitude |E| and (b) longitudinal optical force component F z along the line 
passing through the point with coordinates x — 30 nm and y = 5 nm and parallel to the nanorods for the elliptic (A = 450 nm) 
and hyperbolic (A = 600 nm) regimes of the metamaterial. The insets in panel (a) show the distribution of the electric field 
magnitude in the ££-cut-plane passing through y = 5 nm for the elliptic (A = 450 nm) and hyperbolic (A = 600 nm) regimes 
of the metamaterial. Shaded areas in panel (b) show the region where a pulling force emerges. Electric field amplitude of the 
incident wave is 1 V/m. 


tion are related to the higher multipole contribution and 
the interaction between the particle and the nanorods. 

Lateral force distribution analysis can be provided with 
the following semi-analytical approach. First, the total 
electric field distribution in the nanorod array under ex¬ 
ternal incident wave without the particle is calculated 
numerically with the periodic boundary conditions ap¬ 
plied. Results of the simulation are shown in Fig. [3|d). 
The knowledge of the spatial distribution of the electric 
field magnitude enables calculation of optical forces, with 
two assumptions: (i) the nanoparticle is represented by 
a structureless point electric dipole with a moment // (ii) 
the dipole is assumed to act as a small perturbation to the 
fields of the standalone metamaterial. This means, that 
only collective scattering properties of the nanorod array 
were taken into account, while the mutual re-scattering 
of the field between the particle and nanorods was ne¬ 
glected. Comparison between Figs. |3jb) and (d) verifies 
this approximation both the structure and values of the 
field magnitude are similar. 

The time averaged optical force acting on the point 
dipole is given byP 


<F> = ^V|E| 2 + ^Re[ExH‘], (2) 

where a = a' + ia" is the complex particle’s polarizabil¬ 
ity. The polarizability of the spherical particle is given 


byP 

a = 47reo R 3 — —1 . (3) 

£Au + 2 

The resulting optical force map, calculated using the 
dipolar approximation [Eq. ([2|] appears on Fig. |3jc). 
The arrows show the direction of the force at the cor¬ 
responding point. The color pattern corresponds to the 
absolute value of the force. The remarkable similarities 
between Maxwell’s stress tensor calculations [Fig. !*)] 
and the approximate analytical model [Fig. |3|[c)] sug¬ 
gest the validity of the dipolar model and highlights the 
impact of near-fields on the optical force. It should be 
noted, however, that overall values of optical forces, cal¬ 
culated within those approaches, have about 20% differ¬ 
ence, which is related to the finite size of the particle and 
the mutual field re-scattering between the particle and 
nanorods. 


C. Vertical force component 

The distribution of the lateral optical force component 
(perpendicular to the nanorods) was analyzed in the pre¬ 
vious section. Longitudinal force component F z (parallel 
to the nanorods) is analyzed here. 

As it was already mentioned, the homogenization pro¬ 
cedure averages the near-fields over the unit cell, hence, 
it is inapplicable for estimation of gradient optical force 
in the lateral plane. Nevertheless, the field distribution 
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FIG. 5. (color online) Distribution of the normalized optical force acting on the nanoparticle, situated at the lateral plane 
above the metamaterial. Distance between the center of the nanoparticle and top faces of the nanorods is 10 nm. Arrows show 
the direction of the optical force. Color map shows the distribution of the optical force component parallel to the nanorods 
(F z ). Black solid lines show the the geometrical edges of the nanorods. 


along the z -axis can be roughly estimated considering 
the slab of the nanorod metamaterial as a Fabry-Perot 
resonator in z- direction.^ Therefore, it is reasonable to 
expect a standing wave in the slab (along the z-axis) and 
maxima of the electric field resulting in in-plain trapping 
of the particle. The results of the numerical calculation 
suggest the validity of this hypothesis. Estimation of the 
field maxima position deeply inside the slab can b e pro¬ 
vided by the effective medium approximatiorPlMl but a 
more detailed analysis, that takes into account boundary 
effects, demands numerical simulation. 

The profile of the total electric field magnitude along 
the line parallel to the rods and passing through the point 
x = 30 nm and y = 5 nm calculated without nanoparticle 
for the elliptic (A = 450 nm) and hyperbolic (A = 600 nm) 
regimes is shown in Fig. EFa). The insets show distri¬ 
bution of the total electric field magnitude in xz-pl&ne 
passing through y = 5 nm. One can see that electric 
field distribution along the z-axis strongly depends on 
the wavelength of the incident wave. In the hyperbolic 
regime (A = 600 nm), three distinct field maxima are 
observed - one inside the slab and two in the vicinity 
of its boundaries. In the elliptic regime (A = 450 nm), 
field decays inside the metamaterial and weak oscilla¬ 
tions do not possess sharp field maxima. Additional con¬ 
tribution to those differences (apart from the interplay 
of dispersion regime and geometry, namely Fabry-Perot 
conditions) comes from a strong wavelength dependence 
of losses in goldP^ 

Re(g A u) 

Im(sAu) 


Optical losses cause the reduction in quality factors of 
the modes, smearing out the sharp peaks, as could be 
seen in the case of elliptic dispersion. 

Distributions of the ^-component of the optical force 
along the nanorod for elliptic (A = 450 nm) and hyper¬ 
bolic (A = 600 nm) regimes are shown in Fig. Qb). 
Positions of the stable trapping in transverse planes 
are marked with arrows (note, that the force deriva¬ 
tive should be negative in order to obtain a stable 
equilibrium). Shaded areas on the figure show the re¬ 
gions within the metamaterial, where the optical force 
component F z is directed towards the light source (for 
A = 600 nm). Optical pulling forces or optical attrac¬ 
tion gained considerable attention over the last decade, 
as it provides additional flexible degree of freedom in op¬ 
tical manipulation.^ 

The hyperbolic regime supports three regions of opti¬ 
cal attraction, while the elliptic has only one, as could be 
seen in Fig. 0b). This occurrence could be understood 
as follows: in the elliptic regime both weak gradient of 
the electric field magnitude [see Fig. [4](a)] and high ma¬ 
terial losses of the particle result in the domination of 
radiation pressure [second term in Eq. <@] over the gra¬ 
dient force. The radiation pressure is co-directional with 
the Poynting vector of the incident radiation, so the opti¬ 
cal attraction cannot be obtained (unless negative-index 
materials are involved, which is not the case here). Nev¬ 
ertheless, the first term of Eq. overcomes the second 
one in the vicinity of the nanorod’s edge where strong 
gradient of the electric field magnitude is observed [see 
FigQa)]. In the hyperbolic regime, on the other hand, 
there are several regions where the optical force compo- 


0.3 for A = 450 nm; 
6.2 for A = 600 nm. 
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nent is directed to the light source - that’s the result of 
high quality factor Fabry-Perot modes and dominating 
real part of the particle’s polarizability. 

As a separate case, the particle situated over the meta¬ 
material slab will be considered next. This scenario de¬ 
scribes the case where the metamaterial is used as a sub¬ 
strate for advanced optical manipulation. Results of nu¬ 
merical studies appear in Fig. [5] showing the distribution 
of the vertical optical force acting on the nanoparticle 
in the lateral plane of z = 10 nm above the nanorods. 
It could be seen, that the maximal attraction force on 
the particle emerges in the vicinity of nanorods edges 
(the sample is illuminated from below - see Fig. [5|. Ar¬ 
rows indicate the direction of the optical force. Color 
pattern shows the distribution of the optical force com¬ 
ponent parallel to the nanorods (F z ). Solid white lines 
correspond to F z = 0. Remarkable behaviour of forces 
above the metamaterial substrate could suggest the later 
as an auxiliary nanostructure or metasurface, providing 
additional flexibility on optical manipulation. 

IV. CONCLUSION 

In this work, comprehensive analysis of the optical 
forces acting on a metal nanoparticle placed inside or in 
the vicinity of three-dimensional nanorod metamaterial 
slab was performed. Numerical simulations of finite size 
square unit cells with periodical Floquet boundary con¬ 
ditions enable to take into account all collective effects 
in the metamaterial and estimate optical forces on small 
particles. Unit cells containing 4, 9, and 16 nanorods 
were analyzed and the convergence of the optical forces 
for different positions of the particle was checked. It 
was shown that the smallest unit cell already reproduces 
the effect of optical forces on a particle, situated within 
the infinite metamaterial. Therefore, only four neighbor¬ 
ing nanorods nearest to the particle make the dominant 
contribution to the optical forces. This statement has 


been confirmed with the developed semi-analytical model 
which neglects the particle’s interior and the re-scattering 
effects between the particle and nanorods. Furthermore, 
it was shown that the ’topological transition’ from the 
elliptic to hyperbolic dispersion regime of the metama¬ 
terial, usually having an impact on various light-matter 
interaction processes, is less important for optical forces. 

In-plane optical trapping and optical pulling forces 
were observed. The comprehensive numerical modeling 
enables estimation of optical forces values, normalized to 
incident power and particle’s volume. Values as high as 
2.3x10 s pN/W/nm for both lateral and optical pulling 
forces were predicted. Those results overcome routinely 
reported values (see, e.g., Ref. [36) by an order of mag¬ 
nitude and are comparable with other advanced layouts 
(see, e.g., Ref. EE3- 

The remarkable structure of predicted optomechanical 
interactions (in particular pulling forces), mediated by 
the metamaterial, make the later to be a promising plat¬ 
form for large span on multidisciplinary applications, 
involving demands for precise nanoscale mechanical 
manipulation, including trapping sorting, mixing and 
more. 
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